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ABSTRACT 


Acoustic  propagation  in  the  littoral  regions  of  the  world,  even  over  short  ranges, 
ean  be  eomplex  at  high  frequeneies  (>  5  kHz),  and  applieations  sueh  as  underwater 
deteetion  and  eommunieations  suffer  as  a  result.  To  this  end,  the  Asian  Seas  International 
Aeousties  Experiment  (ASIAEX)  was  condueted  with  funding  from  the  Offiee  of  Naval 
Researeh.  One  phase  of  this  experiment  took  plaee  from  29  May  to  9  June  2001  and 
foeused  on  short-range,  shallow  water  aeoustie  propagation  in  the  East  China  Sea.  This 
thesis  will  be  based  on  some  of  the  measurements  taken  during  the  East  China  Sea 
experiment  by  a  team  from  the  Applied  Physies  Eaboratory  of  the  University  of 
Washington.  The  environmental  parameters  and  array  geometry  used  in  this  and  previous 
modeling  work  are  as  close  as  possible  to  the  experiment. 

The  objeetive  of  this  thesis  is  to  better  understand  the  nature  of  turbulent 
perturbations  and  how  they  affeet  short-range  aeoustie  propagation  in  a  shallow  water 
environment.  This  will  inelude  variations  in  total  energy  of  the  turbulenee,  variations  in 
assoeiated  length  seales,  and  variations  in  depth.  In  addition,  the  influenee  of  other 
variable  faetors  on  signal  eoherenee  will  also  be  examined.  Most  notably,  variations  in 
background  sound-speed  structure  and  souree  depth  variability  will  be  ineluded. 
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I.  INTRODUCTION 


Having  the  ability  to  perform  operations  in  the  littoral  regions  of  the  world  is  one 
of  the  most  erueial  issues  if  a  navy  intends  to  be  fully  eapable.  Sinee  these  regions  differ 
from  the  so-ealled  “blue  waters”  in  many  ways,  even  the  most  thorough  deep  water 
researeh  may  not  be  applieable.  One  of  the  most  extensive  researehes  condueted  in 
littoral  regions  was  the  Asian  Seas  International  Aeousties  Experiment  (ASIAEX), 
funded  by  the  Offiee  of  Naval  Researeh.  The  seeond  phase  of  this  projeet  took  plaee  in 
the  East  China  Sea  (ECS). 

One  eomponent  of  the  ECS  experiment  involved  the  analysis  of  aeoustie  signal 
eoherenee  over  the  length  of  a  vertieal  array  at  relatively  short  range  (500  m)  from  a 
eontrolled  souree.  Eor  eertain  souree/reeeiver  geometries,  direet  path  arrivals,  single 
surfaee  bounee,  and  single  bottom  bounee  paths  eould  be  distinguished,  allowing  analysis 
of  eoherenee  for  those  paths  independently. 

Previous  thesis  work  done  by  ET  Stephen  C.  Karpi  at  the  Naval  Postgraduate 
Sehool  eontributed  to  the  analysis  of  the  direet  path  (no  boundary  refleetions)  arrivals  by 
numerieally  evaluating  the  influenee  of  various  shallow  water  sound  speed  perturbations. 
Karpi’s  researeh  employed  vertieal  eoherenee  to  measure  the  effeets  of  longer  seale, 
internal  wave-like  struetures  and  smaller  seale,  turbulent-like  fluetuations  on  the  direet 
path  propagation.  His  work  showed  that  perturbations  by  an  internal  wave-like  strueture 
affeeted  the  eoherenee  only  slightly.  However,  turbulent-like  fluetuations  eaused 
signifieant  signal  deeorrelation  RMS. 

This  thesis  examines  both  the  model  geometry,  ineluding  souree  and  reeeiver 
depths,  and  the  environmental  variables,  whieh  are  baekground  sound-speed  profde, 
turbulenee  strength  and  the  turbulenee  outer  length  seale.  The  general  geometry  and 
souree  eharaeteristies  used  during  the  East  China  Sea  experiment  by  Peter  Dahl  of  the 
Applied  Physies  Eaboratory  of  the  University  of  Washington  are  the  basis  of  numerieal 
studies  in  this  thesis. 
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The  model  used  in  this  work  is  the  Monterey-Miami  Parabolie  Equation  (MMPE) 
propagation  model.  The  MMPE  model  was  developed  by  Smith  and  Tappert  in  1996.  The 
speeific  version  used  in  this  work  incorporates  the  range-dependent  turbulent  fluctuations 
and  the  source  spectrum  characteristics  employed  in  the  ASIAEX  experiment.  This 
model  provides  an  efficient  algorithm  for  computing  the  broadband  source  response  in 
such  range-dependent  environments. 

A  realistic  turbulence  spectrum  has  multiple  parameters  which  depend  on  various 
features  of  the  ocean  environment.  In  this  work,  only  two  parameters,  a  strength 
parameter  to  adjust  the  overall  RMS  value,  and  an  outer  length  scale  beyond  which 
fluctuations  have  no  contribution,  will  be  considered  variables.  As  a  reference,  10  m  is 
used  for  the  outer  length  scale,  and  four  other  different  values  are  tested.  Because  only 
weak  scattering  is  expected,  the  outer  length  scale  is  not  expected  to  directly  affect  the 
propagation.  However,  the  dominant  contribution  to  the  scattering  is  expected  to  be  at  the 
Eresnel  radius,  and  this  is  on  the  same  order  as  the  reference  scale  of  10  m.  Thus,  we  may 
expect  to  see  some  differences  with  changes  in  the  outer  scale  due  to  the  energy  in  the 
turbulent  perturbation  near  the  Eresnel  scale.  Conversely,  the  turbulence  strength 
parameter  is  expected  to  have  a  direct  effect  on  scattering,  and  thus  affect  the  signal 
coherence. 

Besides  the  two  aforementioned  parameters  of  turbulence,  the  effects  of  variable 
background  sound-speed  profiles  are  examined  in  this  research.  The  reference  sound- 
speed  profile  (SSP)  is  taken  as  the  average  value  of  the  measurements  from  ASIAEX. 
Three  more  SSPs  are  also  examined.  One  of  them  is  very  close  to  the  reference  case,  but 
has  more  variability  along  the  depth.  The  other  two  are  very  different  from  the  reference 
case.  One  of  them  has  a  much  deeper  thermocline  and  steeper  gradient  below,  while  the 
other  has  a  shallower  thermocline  and  generally  less  steep  gradient  below.  The  effects  of 
SSP  are  examined  using  not  only  the  vertical  correlation  but  also  using  some  simple  ray 
tracing  analysis  involving  sound-speed  relative  curvature  (Duda  et  al,  1994). 

The  last  variables  investigated  in  this  work  are  source  and  receiver  depths.  Since 
the  model  computes  the  propagation  at  all  depths,  no  additional  computational  effort  is 
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required  to  obtain  the  solution  at  any  reeeiver  depth.  To  eompute  vertieal  eoherenee, 
short  sub-arrays  are  modeled  by  extraeting  the  solution  over  5  m  vertieal  spans.  The 
reeeiver  sub-array  depths  that  are  analyzed  are  35-40,  45-50,  55-60,  and  65-70  m.  The 
referenee  value  used  for  souree  depth  is  50  m,  while  30,  40,  60,  and  70  m  are  also  tested. 
It  should  be  noted  that  above  30  m  or  below  70  m  of  the  souree  or  reeeiver  depths,  the 
direet  path  propagation  will  overlap  with  those  arrivals  from  surfaee  seattering  and 
bottom  refleetions,  for  the  2  msee  pulse  modeled.  Sinee  the  foeus  of  this  work  is  to  deal 
with  the  direet  path  propagation,  these  depths  are  disregarded  in  the  analysis  range. 

Chapter  II  eontains  an  overview  of  the  theory  behind  the  MMPE  model  and  the 
turbulent  speetrum.  The  implementation  of  these  models  and  the  post-proeessing  methods 
are  also  examined  in  this  ehapter.  Different  model  parameters  and  their  range  of  analysis 
are  eovered  in  Chapter  III.  Chapter  IV  foeuses  on  post-proeessing  and  the  analysis  of 
modeled  data.  The  last  ehapter  presents  eonelusions  derived  from  the  analysis  and 
suggestions  for  future  researeh. 
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II.  NUMERICAL  METHODS  AND  IMPLEMENTATION 


A,  MONTEREY-MIAMI  PARABOLIC  EQUATION  (MMPE)  MODEL 


The  parabolic  equation  (PE)  method  was  introduced  into  underwater  acoustics  in 
the  early  1970’s  by  Tappert.  Due  to  its  efficient  numerical  algorithm  for  computing 
solutions  to  general  range-dependent  problems,  it  has  become  the  most  popular  method 
for  computing  the  acoustic  field  in  complex  underwater  acoustic  environments.  The 
MMPE  model,  developed  by  Smith  and  Tappert  in  the  early  1990’s,  is  based  upon  the 
parabolic  equation  method.  In  this  chapter,  a  simple  derivation  of  the  MMPE  model  is 
described  by  borrowing  heavily  from  Smith’s  descriptions. 

We  begin  by  defining  the  time-harmonic  acoustic  field  in  a  cylindrical  coordinate 
system  by 

P{r,z,(p,(ot)  =  p{r,z,(p)e^'“’‘  .  (1) 


Substituting  Eq.  (1)  into  the  linearized  acoustic  wave  equation  in  cylindrical  coordinates 
leads  to  the  Helmholtz  equation. 


lA 

r  dr 


dp\  1  d^p  d^p 


dr 


+ 


r  d(p  dz 


+  ^^^  +  kln^{r,z,(p)p  =  -A7i;PQd{x-Xs), 


(2) 


where 


=  col Cq  , 

is  the  reference  wavenumber,  and 


(3) 


n{r,z,(p)  =  (4) 

c{r,z,(p) 

is  the  acoustic  index  of  refraction.  In  Eq.  (4),  is  the  reference  sound-speed  and 
c{r,z,(p)  is  the  acoustic  sound-speed.  Note  that  c{r,z,(p)  carries  all  the  information 
about  the  environment,  except  the  density,  which  is  unimportant  in  this  development 
since  we  are  only  interested  in  propagation  in  the  water  column.  Eq.  (2)  implies  that  the 
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The  last  term  in  Eq.  (6)  deereases  quickly  due  to  1/r^  dependence  and  is  generally 
neglected.  The  second  term  accounts  for  azimuthal  coupling  between  different  radials 
and  is  small  enough  to  consider  as  negligible.  Neglecting  this  term  is  referred  to  as  the 
“uncoupled  azimuth  (UNCA)  approximation.” 

The  remaining  Helmholtz  equation  can  be  further  simplified  by  introducing  the 
operator  notations 

and  Q^^={^  +  £  +  iy  ,  (7) 

where 


2  .  15^ 


Eq.  (6)  can  then  be  represented  in  the  form 


(F  +Fap>=o. 

Eactoring  Eq.  (9),  outward  and  inward  fields  are  represented  in  the  form, 
(Pop  +  ^KQop )(Pop  -  iKQop  >  +  iK {Pop .  Qop  ]w  =  0 


(8) 

(9) 

(10) 


Assuming  weak  range  dependence,  the  commutator  \_Pgp,Qgp\  is  negligible. 
Eurthermore,  defining  the  outward  propagating  field  as 


(11) 
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the  outgoing  solution  satisfies 

PVT/ 

Po,^  =  iKQo,^  or  =  (12) 

If  the  baekseattered  energy  is  negligible,  Eq.  (12)  represents  the  eomplete  deseription  of 
the  forward  propagating  aeoustie  energy  in  the  waveguide.  This  equation  is  the 
foundation  for  all  underwater  aeoustie  PE  models. 

To  solve  the  parabolie  equation  given  in  Eq.  (12),  first  the  field  should  be 
deeomposed  into  two  parts,  an  envelope  funetion  and  a  phase  term.  A  reasonable 
deeomposition  for  the  PE  field  funetion  is  defined  by 

‘P  =  .  (13) 

Using  this  modulating  funetion,  we  get  the  PE  for  the  field  funetion  in  the  form 

■  (14) 

PE  solutions  are  based  on  developing  an  approximation  to  the  pseudo-differential 
operator  .  There  are  basieally  three  methods  to  solve  the  PE:  the  Einite  Element  (EE) 

method;  the  Einite  Difference  (ED)  method;  and  the  Split-Step  Eourier  (SSE)  method.  An 
approximation  well  suited  for  the  SSE  algorithm,  on  which  the  MMPE  model  is  based, 
was  introduced  by  Thompson  and  Chapman  (1983)  and  is  commonly  referred  to  as  the 
wide-angle  PE  (T-C  WAPE).  This  approximation  generally  provides  good  accuracy  in 
phase. 

To  obtain  the  SSE  algorithm,  let  us  rewrite  the  PE  in  terms  of  the  Hamiltonian- 
like  operator  as 

,  (15) 

dr 

where 

•  (16) 
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This  Hamiltonian-like  operator  defines  the  evolution  of  the  PE  field  funetion  in  range. 
Solutions  of  the  'F  funetion  at  different  ranges  ean  be  defined  by  a  marehing  algorithm 
of  the  form 


^(r  +  Ar)  =  0(r)^(r)  ,  (17) 

where  0(r)  is  a  unitary  operator  that  progresses  the  solution  out  in  range.  The  propagator 
0(r)  ean  be  approximated  as 

0(r)  «  ,  (18) 

where 


r-fAr 

HoAr)  =  —  \  dr'H^^ir')  .  (19) 

r 

The  advantage  of  the  SSF  algorithm  is  based  on  the  faet  that  is  separated  into 
a  combination  of  spatial  and  differential  operators,  which  can  be  expressed  as 

(20) 

where 


and 


Uo,=- 


\j\  +  s  -  \ 


=  -{\-n) 


(21) 


1  + 


kl  dz^ 


(22) 


Keeping  in  mind  the  assumptions  made  thus  far,  the  propagator  can  be  made  2 
order  accurate  in  range  step  (Ar)  by  employing  the  form. 


nd 


0(r)  = 


Ar 

Ar 

e  ^ 

-i 

e 

_ 

_ 

_ 

(23) 


The  operator  Uop  is  a  multiplication  operator  in  the  z-space  and,  hence,  is  a  diagonal 
matrix.  The  operator  Top  is  not  diagonal  in  the  z-space,  but  it  is  in  wavenumber  space; 

therefore,  the  kz-domain  operator  ^op  is  desired. 
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It  is  now  possible  to  outline  how  the  PE/SSF  algorithm  steps  out  in  range.  To 
begin,  the  PE  field  funetion  is  defined  at  some  range,  r,  in  the  z-domain.  This  funetion  is 
then  multiplied  by  the  operator,  e  ^  defined  at  the  beginning  of  the  range  step. 

A  transformation  to  the  ^^-domain  then  oeeurs  so  that  we  ean  multiply  by  e  .  This 
result  is  then  transferred  baek  to  the  z-domain  to  be  multiplied  by  the  z-spaee  operator 
defined  at  the  end  of  the  range  step,  e  diserete  fast  Fourier  Transform 

(EFT)  is  being  used  in  the  eode  assuming  the  convention 


y,{z)  =  FFT{y>{k^)) 


(24) 


and 


y>{K)  =  IFFT{y,{z))  . 
The  PE/SSF  implementation  can  then  be  represented  by 


(25) 


Ar,z)  =  e 


Af 

-ik^ — (r+Ar,z) 

”2  VfT 


(  i.  i  ^ 

-tK-rUop(r,^)  ,  , 

e  y/[r,z) 

^-ik^^Arf^p  {k^ )  ^  JppJ' 

V  )_ 

(26) 


where 


Uoj 


(27) 


B,  TURBULENT  SOUND-SPEED  PERTURBATION  THEORY 

Acoustic  propagation  through  a  medium  is  subject  to  the  properties  of  that 
medium.  When  a  spherical  wavefront  from  a  point  source  traverses  a  medium  containing 
fluctuations  in  sound  speed,  portions  of  the  front  travel  faster  or  slower  than  other  parts 
of  the  wavefront  and  the  front  distorts  as  a  result.  In  the  ocean,  such  sound-speed 
inhomogeneities  result  from  phenomena  like  layering,  internal  waves,  and  turbulence 
(Duda,  et  al,  1988).  In  this  work,  it  is  assumed  that  perturbations  of  sound  speed  are 
generated  purely  by  turbulent  fluctuations.  These  fluctuations  are  characterized  using 
statistical  relationships  developed  in  previous  work  (Tatarskii,  1961;  Duda,  et  al,  1988). 
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The  statistical  theory  of  turbulence  characterizes  the  turbulent  perturbation  field  as 
realizations  of  a  random  function.  Such  a  characterization  is  only  valid  when  changes  in 
signal  amplitude,  due  to  fluctuations  in  the  medium,  do  not  exceed  the  average  amplitude 
of  the  signal.  This  is  referred  to  as  the  “weak  scattering  regime”  (Duda  et  al,  1988). 

If  we  define  the  index  of  refraction  perturbation  as 

Sn{r)  =  ^^^  =  ^{r),  (28) 

c{r) 

then  the  variance  of  /u{r)  is  related  to  the  spectral  density  of  the  perturbations,  S ^  (K) , 
by 


var(//(F))  = 


1 


\s^iK)d^K, 


(29) 


where  the  integral  is  a  three-dimensional  integral  over  all  wavenumber  space.  The 
variance  is  defined  by 


so 


var(//(F))  =  (|//(F)n  =  — - I  u\f)d^r  ,  (30) 

'  '  volume 


volume 


\s^{K)d^K. 


If  we  define  the  spectrum  (transform)  of  ju{r  )  as  FJ^JC)  by 


//(r)  = 


f  F(K)/^*^d^K, 

{Inf  I 


then  the  correlation  of  jLi{r)  satisfies 


(31) 


(32) 


00  1  °°  _  2  - 

J  MF)/(F-FVV==— y  J  d^K\F^{K)\  .  (33) 

_O0  ^  /  _00 

I  “  1^ 

Thus,  \F^{K)\  is  the  Fourier  transform  of  the  correlation  function.  It  is  easy  to  see  that 
by  setting  r '  =  0  ,  the  square  of  the  signal  spectrum  and  the  spectral  density  are  related  by 


10 


It  follows  then  that 


=  {volume)* S^{K)  . 


(34) 


\F^{K)\  =  \^olume*  S ^{K)~^  .  (35) 

Note  that  //(F) is  unitless  while  F^{K)  and  (^) both  have  units  of  m^. 

Similar  reasoning  holds  for  a  two-dimensional  shoe  at  a  fixed  value  of  one 
variable,  e.g.  n{x,y  =  0,z) ,  sueh  that 

-«  <  00  00 

- - JJ 7  =  0,  z)dxdz  =  J  J  {K^ ,  )dKJK^ ,  (36) 

where  {K^ =  {k^ ,k^,ky  =0)  is  the  two-dimensional  speetral  density  along  the 

slieey  =  0.  Furthermore, 

00  00 

Mix,  z)  =  — ^  J  J  i^x ,  K ,  (37) 

V'^  ^  7  —00  —00 

where  G^{k^,k^)is  the  speetrum  (transform)  of  the  two-dimensional  perturbation  field.  It 
follows  that 

\G/k,.k,)\^[Area»r^{k,.k,)]"\  (38) 

The  two-dimensional  speetral  density  ean  be  obtained  from  the  three-dimensional 
speetral  density  aeeording  to  (Tatarskii,  1961), 

1 

r,(KX)  =  —\s^(K,k^,k,)dk,.  (39) 

-00 

At  this  point  in  the  turbulent  analysis,  perturbations  in  sound  speed  have  been 
related  to  the  speetral  density  of  the  turbulent  fiuetuations.  Developing  an  expression  for 
the  turbulent  perturbation  field  and  ealeulating  its  statisties  remain  to  be  done.  Tatarskii 
approximates  the  turbulent  field  as  isotropie  and  homogeneous  by  using  the  method  of 
strueture  funetions  (Tatarskii,  1961).  Aeoepting  that  large-seale  inhomegeneities  may  be 
present,  strueture  funetions  restore  homogeneity  by  redueing  the  seale  of  the  analysis. 
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Standard  statistical  relations  are  then  applied,  resulting  in  a  three-dimensional  spectral 
density  function  of  the  form  (Tatarskii,  1961), 

E{K)  =  A  ^  ,  (40) 

where  ,  A  is  a  sealar  multiplier,  and  e  is  the  energy  dissipation  rate. 


In  the  literature  (e.g.,  Duda,  et  al,  1988),  this  speetral  density  has  been  adapted  as 
follows.  To  avoid  the  speetral  density  beeoming  infinite  at  if  =  0,  a  wavenumber 
threshold  (outer  scale)  limit  is  defined  as  kt  such  that 


EiK)^ 


1 


(K^+k^) 


2\ll/6 


(41) 


thus  maintaining  the  K  power  law  from  Tatarskii’s  speetrum.  A  highpass  filter  of 
the  form 


is  also  added  to  force  the  spectral  density  to  zero  at  K=0. 


(42) 


For  high  wavenumbers,  a  high-frequeney  eut-off  related  to  the  Batchelor 
spectrum  is  given  by 


where 


Rb(K)  =  qxp 


(43) 


A  ^  A 


1/4 


kg  — 


j 


(44) 


is  the  Batchelor  wavenumber,  e  is  the  depth  averaged  kinetic  energy  dissipation  rate  in 
W/kg,  V  =  1.40e-6  m  /s  is  the  kinematie  viscosity,  and  Kj  ~1.00e-7  m  /s  is  the  thermal 
diffusivity  of  sea  water.  The  eonstant  q  is  an  order  unity  faetor,  and  3.7  is  used  here  in 
concordance  with  earlier  work  (Duda,  et  al,  1988). 
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The  depth  averaged  kinetic  energy  dissipation  is  defined  as  (Henyey  et  al,  1986), 


e=0.3  — ( 


[  4  ^ 

2 

cosh  ' 

r  Niz)^ 

\L^  ' 

y  fi  J 

(45) 


The  brackets  (  indicate  a  depth  average,  j»  =  1  is  the  characteristic  mode  for  shallow 
water,  bE^j^  =  0.5m  is  a  measure  of  internal  wave  intensity,  fj  is  the  inertial  frequency 
in  rad/sec,  and  the  buoyancy  frequency  N{z)  may  be  computed  from 


N{z)  =  g^'^ 


1 


dz 


1/2 


(rad/  sec). 


(46) 


where  p{z)  and  (z)  are  the  density  and  potential  density,  respectively,  as  a  function 

of  depth.  This  expression  for  "e  is  valid  for  turbulence  created  by  internal  wave  decay 
(Henyey  et  al,  1986).  It  is  an  expression  developed  for  deep  water  and  ignores  other 
possible  sources  of  turbulence. 


By  combining  these  quantities,  the  three-dimensional  spectral  density  has  the 

form 


d) 

S(K)  =  — ,  - v/?«(^)  • 

{K^ +kf  "" 


(47) 


From  this  expression,  the  two-dimensional  spectral  density  is  calculated  below.  The 
integration  is  applied  only  to  the  portion  defining  the  turbulent  spectrum  (the  scaling 
factor  and  the  wavenumber  roll-off  terms  are  added  afterwards). 


®3Z2(^)  (^2^^2yi/6- 


(48) 


The  two-dimensional  spectral  density  is  then 


^  CO 

^2Dik,,k^)  =  —  \Q>,^{k^,k^,k^)dky 


1 

J-f 

“)  TT  J 


I 


27riik:+k  +k:+k:y 


rdk,. 
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(49) 


1  V^r(i) 

27r[kl+kl+k^Y"  m) 

Employing  the  same  highpass  filter,  high-wavenumber  Batchelor  spectrum  cutoff, 
and  turbulent  strength  parameter  scaling,  the  two-dimensional  spectral  density  of 
turbulent  fluctuations  of  the  index  of  refraction  is 


VAK)  = 


r(l)  4 


In  r(^)  {K^+k^y^^  K^+kf 


R,(K). 


Note  that  V  (K)  has  units  of  m^  as  expected. 


(50) 


C.  IMPLEMENTATION  METHODS 

1.  Numerical  Implementation  of  Turbulence 

In  the  above  discussion,  V^(K) ,  the  two-dimensional  spectral  density  of 
refractive  index  perturbations,  was  related  to  /j(x,  z) ,  the  perturbation  field.  A  random 
realization  of  the  two-dimensional  turbulent  field  of  sound-speed  perturbations  is 
obtained  by  employing  the  two-dimensional  spectrum,  {k^  ,ky  ,  and  defining 

..  00  00 

A(az)  =  -^  I  I  [g  {k^,kyB{k^,kye^^’‘^’''^^^^  (51) 

(2;z-)  J 

where  B{k^,k^)  and  6{k^,k^)  are  two-dimensional  random  number  functions.  Since 

the  complex  amplitude  of  each  wavenumber  component,  Be^ ,  should  exhibit  a  normal 
distribution,  the  random  phase  of  each  component  is  a  uniformly  distributed  random 
variable  on  the  unit  circle.  In  other  words,  we  define 

6  =  2nr^ ,  (52) 

where  is  a  uniformly  distributed  random  variable  in  the  interval  [0,1].  The  magnitude 
B,  however,  exhibits  a  Rayleigh  distribution.  Because  we  are  in  fact  generating  a 
realization  of  the  spectral  density,  we  must  consider  the  magnitude-squared,  which  has  a 
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negative  exponential  distribution.  Thus,  the  random  amplitude  of  each  wavenumber 
component  is  obtained  by 


B  =  (53) 

where  Tj  is  another  independent,  uniformly  distributed  random  variable  in  the  interval 

[0,1]. 


Since  we  are  implementing  this  in  a  model  with  a  discretely  sampled  grid  of  the 
environment,  we  must  consider  discrete  Fourier  transforms  computed  using  an  FFT 
algorithm.  Thus,  with  scaling  appropriate  to  the  FFT  algorithm  being  employed. 


(27rf 


M  N 

Z  Z  [g,  (*„ .  K.  )5(*,, .  K.  ’  ] 

m=l  n=\ 


e 


(54) 


From  the  previous  analysis  of  the  relationship  between  the  perturbation  spectrum  and  the 
spectral  density. 


where  we  here  define  the  spectral  density  up  to  a  scaling  factor 


ZAK)  = 


r(l)  i 


R,(K). 


-2 


(56) 


Itt  r(^)  (K^+k^y'^  K^+kf 

Note  that  is  the  fixed  turbulent  scaling  parameter  and  may  be  adjusted  to  produce  the 
desired  RMS  perturbation. 

Combining  the  above,  we  find 


\K  AjL  _  M  N 


i(K, 


-KnH) 


{In) 


n=\  n=\ 


=  y/AreaFFT^°  (K)B(K)e‘''^^^  ]  (57) 

For  a  given  computational  area,  the  factor  Area  =  RZ,  where  the  turbulent  field  exists. 
Furthermore, 
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(58) 


so 


Thus 


^k^^ky 


Area 


Area  = 


Ak^Ak 

A  y 


(59) 


M{X.,Z,)  =  ^Ai^FFT^^  [v^J\K)B{K)e^‘^^‘^^\  (60) 

Once  a  realization  of  //(x,  z)  is  computed,  we  simply  multiply  by  Co  to  obtain  the 
sound-speed  perturbation, 

5c{x,z)  =  c^iu{x,z) .  (61) 

Technically,  we  might  multiply  by  c(x,z),  a  sound-speed  profile  with  range.  However, 
variations  in  c(x,z)  are  typically  less  than  two  percent  and  so  using  Co  is  a  reasonable 
approximation.  Furthermore,  note  that  we  are  employing  the  perturbation  in  eylindrieal 
eoordinates,  replacing  x  by  r,  so 

Sc{x,z)  ^  Sc{r,z) .  (62) 

This  is  also  justified  so  long  as  we  only  eonsider  propagation  along  the  two-dimensional 

slice  in  (r,z).  See  Figure  1  for  an  example  realization  of  the  sound-speed  perturbation 
field. 
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Realization  of  Turbulent  Sound  Speed  Perturbation 


Range  (m) 


Figure  1.  Sample  realization  of  the  turbulent  sound  speed  perturbations  in  m/s.  The  RMS 
sound  speed  fluctuation  is  0.41  m/s. 


To  compute  the  parameters  A^(z)and£  ,  we  need  to  extract  both  the  buoyancy 
frequency,  N(z),  and  the  variance  of  the  temperature  gradient,  ^  The  former 

may  be  defined  as 


N{z)  =  g 


1/2 


1  5/7„(z) 


1/2 


(rad/  sec), 


(63) 


p{z)  dz 

where  p{z)  and  Pp{z)  are  the  depth  profiles  of  the  density  and  potential  density, 

respectively,  and  g  is  the  acceleration  of  gravity.  Both  densities  were  measured  during 
the  ASIAEX  experiment  by  CTD  casts. 
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2. 


Creation  of  Source  Spectrum 


During  the  experimental  periods  of  interest,  charaeteristic  eontinuous  wave  (CW) 
and  frequency  modulated  (FM)  pulses  were  transmitted.  The  CW  pulses  of  2  to  3  msec  in 
duration  had  center  frequencies  of  4,  8,  16,  and  20  kHz.  The  FM  pulses  are  not  analyzed 
here.  To  simulate  a  source  pulse  of  2  msec,  an  analysis  must  be  performed  in  the  time 
domain  after  the  model  has  run.  A  band-pass  filter  of  approximately  8  kHz  must  also  be 
incorporated,  to  be  consistent  with  the  work  being  done  at  APL-UW.  However,  before 
these  concerns  can  be  addressed  in  detail,  a  brief  summary  of  time  domain  processing 
within  the  MMPE  is  necessary. 

Recall  from  the  development  of  the  parabolic  equation  that  the  time-harmonic 
acoustic  field  (ignoring  azimuthal  dependence)  is  given  by 

P{r,z,o)t)  =  (64) 

A  discrete  FFT  algorithm  is  used  to  represent  the  field  in  the  time  domain  as 

P{r,z,t)  =  FFT[p^{r,z)\  =  Y,pSr,z)e“'‘  ■  (65) 

CO 

In  other  words,  Eq.  (64)  represents  a  single  chosen  frequency  and  a  single  component  of 
the  time  dependent  field.  Thus,  for  analysis  in  the  time  domain,  the  model  must  be  run  at 
multiple  frequencies.  The  inputs  required  by  the  model  for  such  analyses  are  a  center 
frequency  ifo),  a  bandwidth  {BW),  and  the  number  of  frequency  bins  to  compute 
(transform  size,  N).  In  this  research,  the  bandwidth  used  was  8184  Hz  with  a  transform 
size  of  1024,  consistent  with  Karpi’s  previous  work. 

To  view  the  arrival  structure  at  the  receive  array  (r  =  R,  where R  is  the  range  of 
the  array  from  the  source),  we  must  compute  the  complex  field  p^  ^iz)  for  many 

frequencies  and  then  Eourier  transform  to  obtain  Pj^{z,t)  ,  the  complex  pressure  values  in 
time/depth  space.  An  example  P,^{z,t)  is  shown  in  Eigure  2.  Outside  the  given  bandwidth 
of 
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(66) 


f  , 

(  _  BW\ 

</< 

( 

+ 

1  2  ; 

^  2  ; 

the  source  amplitude  is  zero.  The  MMPE  model  ignores  the  zeroes  outside  the  bandwidth 
and  computes  a  base-banded  field  to  avoid  unnecessary  computations.  The  resulting 
expression  is 


where  =lnf^.  The  complex  pressure  values  ^ (z)  are  thus  placed  in  frequency 
bins  symmetrically  about  co-co^  =0  with  limits  of  ±^BJV  . 


Transmission  Loss  (dB  re  1m) 


Figure  2.  Pulse  arrival  structure  for  a  source  depth  of  50m,  RMS  sound  speed 
perturbation  of  about  1  m/s,  and  an  outer  cutoff  length  of  10m. 

Note  the  transform  pair  actually  used  is 
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(68) 


icoR/ 


where  is  the  PE  field  function  and  the  overall  phase  factor  e  is 

neglected.  As  a  result,  arrival  times  are  given  as  values  of  “reduced  time”,  r,  such  that 


r  = 


^  r'^ 

t - 

V 


(69) 


where  is  the  reference  sound  speed.  The  complex  pressure  values  P'i^{z,t)  are  thus 


placed  in  time  bins  symmetrically  about  r  =  0,  with  limits  of  ’  where  T  is  the  total 
length  of  the  time  window.  The  time  step  is  given  by 


At  = - «  .122msec  , 

BW 


(70) 


with  a  corresponding  frequency  step  of 


A/  = 


BW 


N-\ 


^  =  8Hz. 
1024-1 


(71) 


Running  the  model  at  multiple  frequencies  has  thus  enabled  the  analysis  of  the  time 
dependent  field. 


The  APT  team  from  the  University  of  Washington  attempted  to  remove  some  of 
the  out-of-band  noise  from  the  ASIAEX  data  by  applying  a  digital  band-pass  filter  to  the 
signal.  During  the  filtering  process,  they  preserved  8  kHz  of  bandwidth,  using  roughly  a 
+/-  4  kHz  band  around  the  center  frequency.  However,  the  band-pass  filter  was  not 
always  applied  symmetrically  about  the  center  frequency  of  the  signal.  In  order  for  the 
model’s  results  to  be  comparable,  the  model  was  run  using  a  center  frequency  defined  by 
the  center  of  the  filter.  For  example,  on  a  model  run  meant  to  simulate  a  2  msec  pulse  at 
20  kHz,  a  center  frequency  of  19  kHz  was  input  to  the  model,  corresponding  to  the  center 
of  the  filter.  Thus,  the  output  data  is  centered  properly,  but  still  representative  of  a 
continuous  source. 


20 


To  simulate  a  2  msec  pulse,  a  filter  function  was  applied  to  the  data.  A  square 
filter  function  (shown  below  in  Figure  3)  was  chosen  to  represent  the  pulse  in  the  time 
domain.  The  length  of  the  filter  corresponds  to  the  number  of  discrete  time  steps  that  are 
needed  to  comprise  the  pulse,  found  by 
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Figure  3.  Square  filter  function  representing  a  2  msec  pulse. 


2msec  ,  , 
= - =  16, 


pulse 


At 


(72) 


where  the  time  step  is  0.122  msec,  as  shown  in  Eq.  (70).  This  square  filter  will  then  be 
transformed  to  a  source  spectrum  in  the  frequency  domain,  as  shown  in  Figure  4,  and 
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applied  to  the  output  data.  However,  the  output  data  is  centered  at  the  middle  of  the 
bandpass  fdter,  and  the  source  spectrum  should  be  applied  at  the  center  frequency  of  the 
source.  Thus,  a  shift  may  be  required.  In  the  above  example,  a  shift  of  1  kHz  would  be 
needed  to  center  the  source  spectrum  at  20  kHz  instead  of  19  kHz.  Table  1  provides  a  list 
of  the  fdters  used  and  the  corresponding  shift  required. 


Experimental 
Carrier  Freq.  (kHz) 

Minimum  Freq.  of 
Filter  (kHz) 

Maximum  Freq.  of 
Filter  (kHz) 

Center  Freq.  of 
Filter  (kHz) 

Required  Shift 
(kHz) 

4 

2 

10 

6 

2 

8 

5.5 

13.5 

9.5 

1.5 

16 

12 

20 

16 

0 

20 

15 

23 

19 

1 

Table  1.  Characteristics  of  the  bandpass  filters  used  in  data  processing  for  each 
experimental  carrier  frequency.  The  required  shift  is  the  difference  between  the  center  of 
the  filter  and  the  experimental  carrier  frequency. 


To  allow  for  shifting,  the  generated  source  spectrum  must  have  a  wider  bandwidth 
(BfV)  than  the  output  data.  The  bandwidth  is  simply  doubled,  resulting  in  a  new  time  step. 
At' ,  given  by 

At' =  — - —  =  — - — »  0.061msec,  (73) 

BJV'  2BW 


and  thus  twice  as  many  steps  are  required  for  a  2  msec  pulse. 


,  2msec 

n,uise  -  -32. 


(74) 


Note  that  the  total  time  window  T  and  the  frequency  step  Af  remain  constant.  The  end 
result  is  an  identical  spectrum  with  more  frequency  information  added  on  each  side  to 
make  the  spectrum  twice  as  long.  The  original  spectrum  can  be  seen  below  in  Figure  4 
and  Figure  6  illustrates  the  extended  version. 
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Source  spectrum 
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Figure  4.  Source  spectrum  created  from  2  msec  square  window  with  a  bandwidth  of  8184 
Hz. 


In  the  frequency  domain,  the  square  window  takes  the  expected  sine  function 
form,  given  by 


sinc(x)  = 


sin(x) 


X 


(75) 


The  source  spectrum  developed  above  must  now  be  combined  with  an  appropriate  band¬ 
pass  filter.  A  rectangular  window  with  a  Hanning  taper  at  each  end  was  chosen  for  this 
purpose.  The  Hanning  window  that  was  used  can  be  defined  as 


f 

6  k  'W 

1-cos 

In 

V 

\  n  +  lyj 

where  k=  l,...,n. 


(76) 
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The  shape  of  a  Hanning  window  can  be  seen  in  Figure  5.  The  Hanning  window 
was  subsequently  divided  into  halves,  each  half  being  applied  to  an  end  of  the  rectangular 
window,  and  the  total  filter  width  being  equal  to  8  kHz  in  this  case  (spanning  15-23  kHz). 
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Figure  5.  Typical  Hanning  Signal-Processing  Window 

This  system  response  filter  was  then  used  as  a  multiplier  on  the  2  msec  pulse 
spectrum,  thus  completing  the  source  amplitude  function.  An  example  bandpass  filter 
and  source  spectrum  are  shown  together  in  Figure  6  for  a  carrier  frequency  of  20  kHz. 


Hanning  Window 


j _ I _ I _ di. 
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Source  Amplitude  Function 


Figure  6.  Extended  source  spectrum  overlaid  with  8  kHz  bandpass  filter.  In  this  case,  a 
carrier  frequency  of  20  kHz  is  used,  thus  requiring  a  1  kHz  shift. 

3,  Direct  Path  Extraction 

At  the  boundaries  of  the  medium,  the  acoustic  signal  interacts  with  the  surface  of 
the  water  and  the  iso-speed  bottom  (set  at  1600m/s).  In  order  to  remove  the  interface 
effects  from  the  analysis  and  focus  strictly  on  the  volume  effects,  the  direct  path  arrival 
was  extracted.  As  can  be  seen  in  Figure  7,  for  a  source  depth  of  50  meters,  several 
propagation  paths  are  evident  in  the  pressure  time  series  (received  at  500m  range). 
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Arrivals  at  Shallow  Array,  Lo  =  10  meters,  rms  =  1.125m/s 


xIO’^ 


Figure  7.  Pressure  time  series  received  at  shallow  and  deep  depths  at  a  range  of  500 
meters.  Lo  is  the  outer  cutoff  length,  and  RMS  is  the  RMS  value  of  the  sound  speed 
perturbation. 

A  source  depth  of  25  meters  in  some  cases  results  in  a  direct  path  signal 
indistinguishable  from  the  surface  bounce  arrival.  As  in  earlier  research  by  Karpi,  a  50- 
meter  source  depth  will  be  considered  as  a  reference.  Plots  similar  to  Figure  7  were 
analyzed  for  each  model  run  to  determine  the  times  encompassing  the  direct  path  arrival. 
The  time  series  was  then  extracted  directly,  without  filtering,  as  shown  in  Figure  8.  The 
direct  path  extraction  is  performed  at  six  varying  depths  across  each  of  the  roughly  5  m 
length  vertical  array  clusters.  The  uppermost  depths  of  each  sub-array  were  defined  as 
the  numerical  mesh  points  closest  to  35,  45,  55,  and  65  m.  The  depths  of  the  five  lower 
elements  were  then  approximately  15,  45,  105,  225,  465  cm  below  the  uppermost 
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element,  thus  spanning  roughly  5  m.  Therefore,  24  total  signals  are  extracted,  one  for 
each  of  the  six  elements  in  the  four  sub-arrays. 


Figure  8.  Extraction  of  the  Direct  Path  Arrival. 

4,  Vertical  Spatial  Correlation 

It  is  now  possible  to  obtain  a  quantitative  estimate  of  signal  degradation  by  cross- 
correlating  the  extracted  direct  path  signals.  Correlation  is  a  measure  of  the  similarity 
between  two  waveforms.  In  this  case,  signals  recorded  at  different  depths  along  the  same 
propagating  wavefront  are  compared.  If  signal  degradation  was  not  present,  such  as  when 
propagating  through  a  homogeneous  medium,  we  would  expect  the  received  pressure 
signal  to  be  very  similar  at  different  points  along  the  front.  As  the  similarity  of  the 
received  signals  along  the  front  decreases  due  to  inhomogeneities,  we  say  that  our 
wavefront  has  degraded. 

In  the  experiment,  each  of  the  six  elements  in  a  vertical  array  cluster  has  a 
corresponding  received  pressure  signal  recorded  at  the  depth  of  the  element.  Thus,  15 
different  cross-correlations  are  possible  (element  1  with  element  2,  element  2  with 
element  3  and  so  forth),  in  addition  to  the  auto-correlation,  each  with  a  depth  separation 
given  by  the  separation  of  the  elements  being  correlated.  The  final  output  is  a  graph  of 
the  normalized  maximum  vertical  cross-correlation  as  a  function  of  the  depth  difference 
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of  the  elements  being  correlated  with  16  points.  These  16  points  include  each  possible 
combination  between  elements  and  the  autocorrelation  of  the  signal  at  element  one. 


Practically,  each  extracted  direct  path  signal  is  cross-correlated  using  the 
MATLAB  function  xcorr,  which  generates  one  number  for  the  normalized  max 
correlation.  Note  that  “normalized”  means  the  autocorrelation  (the  correlation  of  the 
signal  with  itself)  is  assigned  a  value  of  1.  The  normalized  cross-correlation  can  be 
represented  as  follows: 


Cp,,p2(j)  =  e 


\pi  \T)p^\T  -T)dT 

'  A  \T) 

Ut\ 

A  \T) 

"dT 

(77) 


where  and  p^*  are  the  first  base-banded  pressure  signal  and  the  complex  conjugate  of 
the  second  base-banded  pressure  signal,  respectively. 
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III.  MODEL  VARIABLES 


A,  EXPERIMENTAL  AND  MODEL  GEOMETRY 

The  input  parameters  to  the  MMPE  model  were  taken  mostly  from  data  sets  of  the 
East  China  Sea  experiment,  and  the  geometry  and  source  characteristics  were  modeled 
after  the  short-range  experiment  conducted  by  the  group  from  APE-UW.  The  horizontal 
range  between  the  source,  placed  at  a  depth  of  either  25  or  50  m,  and  receivers  was  set  at 
500  m,  and  the  water  column  depth  was  set  at  100  m.  In  the  ASIAEX  experiment,  two 
four-element  vertical  arrays  were  placed  at  nominal  depths  of  26  and  52  m.  The 
distances  between  the  elements  of  each  array  were  13,  30,  and  60  cm.  Eigure  9  below 
shows  the  experimental  configuration. 
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Eigure  9.  Experimental  Geometry  Ifom  ASIAEX 
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The  MMPE  model  is  capable  of  incorporating  both  shallow  and  deep  bottom 
characteristics.  In  this  work,  however,  the  deep  bottom  depth  was  set  at  3000  m,  deep 
enough  not  to  affect  propagation.  The  sound  speed  of  the  shallow  bottom  was  set  at  a 
constant  of  1600  m/s.  This  version  of  the  MMPE  model  requires  a  reference  sound  speed, 
Co,  set  at  1500  m/s  for  this  research.  It  should  also  be  noted  that  all  results  presented  here 
are  for  runs  done  in  “efficiency”  mode  as  opposed  to  “accuracy”  mode.  The  difference 
between  the  two  modes  is  in  mesh  size,  Ar ,  Az  .  Efficiency  mode  was  chosen  because  it 
is  much  faster  than  accuracy  mode,  and  in  most  cases  should  be  accurate  enough  for  our 
direct  path  analysis. 

Though  this  research  is  mostly  concerned  with  the  effects  of  environmental 
variables,  it  is  also  a  point  of  interest  to  learn  about  the  effects  of  the  geometric  variables, 
such  as  source  and  receiver  depths.  Since  the  model  computes  the  solution  at  all 
computational  depth  grid  points,  it  is  a  simple  matter  to  analyze  the  dependence  of  signal 
coherence  on  receiver  depth.  Thus,  the  numerical  work  presented  here  differs  in 
geometry  from  ASIAEX  in  that  a  total  number  of  four  receiver  sub-arrays  with  six 
elements  each  are  processed.  All  analysis  in  this  work  used  uppermost  sub-array  element 
depths  of  35,  45,  55,  and  65  m,  with  total  apertures  of  roughly  5  m. 

1,  Source  Depth 

In  contrast  to  the  ASIAEX  geometry,  the  influence  of  live  different  source  depths 
is  examined  in  this  research.  A  depth  of  50  m  is  assumed  as  the  reference  source  depth. 
Other  source  depths  included  30,  40,  60,  and  70  m.  The  reference  geometry  for  the  50  m 
source  depth  and  the  four  sub-arrays  is  displayed  in  Eigure  10. 


30 


Air/ Water  Interface 


500  meters 


100  meters 


50  meter 
Source 


Each  Cluster 

— > 

consist  of 

6  elements 

15  cm 

35  meter 

,30  cm 

Receiver  | 

1 

60  cm 

Cluster 

N 

120  cm 

1 

1 

N 

240  cm 

Iso-speed  bottom,  c  =  1600  m  /  s 


Not  to  Scale 


Figure  10.  Model  Geometry  Used  in  This  Research 


B.  ENVIRONMENTAL  MODELS 

The  focus  of  this  research  is  on  the  influence  of  turbulence,  so  parameters 
associated  with  the  description  of  the  turbulent  field  will  be  varied  to  determine  their 
influence.  Specifically,  we  will  examine  the  effect  of  variable  outer  length  cut-off  scale, 
kt,  and  overall  turbulent  strength  scaling  parameter,  (j)!.  These  parameters  have  a  direct 
effect  on  the  structure  of  the  turbulence,  and  thus  a  direct  effect  on  the  coherence  of  the 
propagation.  In  addition  to  these  parameters,  we  shall  also  examine  how  source/receiver 
geometry,  which  affects  propagation  path,  may  influence  the  signal  coherence.  Since  the 
propagation  paths  are  related  to  the  background  sound  speed  profile,  analysis  will  also  be 
performed  with  other  SSPs. 
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1.  Turbulence  Strength 

As  noted  in  Chapter  II,  the  spectral  density  of  turbulence,  given  by  Eq.  (65),  has 
an  overall  scaling  factor,  ,  which  we  refer  to  as  the  turbulent  scaling  parameter.  Any 

desired  RMS  perturbation  may  be  obtained  by  adjusting  this  scaling  parameter.  In  this 
work,  the  turbulent  scaling  parameter  has  been  fixed  at  live  specific  values,  considered  to 
range  over  realistic  turbulence  perturbation  scales  from  very  weak  to  very  strong.  One 
value  was  considered  a  reference  value,  two  values  were  examined  below  reference,  Ref- 
1  and  Ref-2,  and  two  values  were  examined  above  reference,  Ref+1  and  Ref+2.  Sample 
values  for  the  RMS  perturbation  for  each  of  the  turbulent  strength  cases  is  provided  in 
Table  2. 


Case 

Ref-2 

Ref-1 

Reference 

Ref+1 

Ref+2 

RMS 

0.112 

0.224 

0.448 

0.897 

1.793 

Table  2.  RMS  Perturbation  of  Sound  Speed  (m/s)  for  Different  Turbulent  Strengths 

2,  Turbulence  Outer  Length  Scale 

The  outer  length  cutoff  scale,  k, ,  in  Eq.  (65)  serves  as  a  low  wavenumber  cutoff, 
limiting  contributions  from  the  turbulent  perturbation  to  wavenumber  values  roughly 
greater  than  k, .  Since 


(78) 


then  \  is  a  wavelength  threshold  that  limits  the  spatial  scale  of  the  turbulent  fluctuations. 


The  weak  scattering  theory  dictates  that  the  scattering  of  a  propagating  signal 
results  mainly  from  perturbations  on  the  scale  of  the  Eresnel  radius.  The  Eresnel  radius, 
Rf ,  is  given  by  (Duda  et  al,  1998) 

,  (79) 
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where  r  is  the  transmission  range  and  A,  is  the  wavelength  of  the  signal.  For  a  signal  with 
a  frequency  of  4  kHz,  the  wavelength  is  given  by 


;l  = 


/ 


1500m/sec 

4000Hz 


0.375  m. 


(80) 


Thus  for  a  transmission  range  of  500  m,  the  Fresnel  radius  is 

^/,5oo«,,i6mz  =  V500m*0.375m  □  13.7m  .  (81) 

The  corresponding  “Fresnel  wavenumber”  is  then  about  0.46  m"'.  Table  3  shows  the 
Fresnel  radius  and  Fresnel  wavenumber  values  for  frequencies  examined  in  this  research. 


Frequency  (kHz) 

Wavelength  (m) 

Fresnel  Radius  (m) 

Fresnel  Wavenumber  (1/m) 

4 

0.375 

13.693 

0.459 

8 

0.188 

9.682 

0.649 

16 

0.094 

6.847 

0.918 

20 

0.075 

6.124 

1.026 

Table  3.  Fresnel  Radius  and  Fresnel  Wavenumber  Values  for  Different  Frequencies  at  a 
Range  of  500  m 

Though  the  values  for  the  Fresnel  radius  vary  between  6  and  14  m,  10  m  is  used 
as  the  reference  case  for  the  turbulence  outer  length  scale  for  all  frequencies  in  the 
research.  The  other  values  examined  in  this  work  are  2  m  and  5  m,  as  Ref-2  and  Ref-l, 
respectively,  and  25  m  and  50  m,  as  Ref+1  and  Ref+2,  respectively.  As  the  scale  value 
was  increased  from  low  values,  the  most  distinguishable  change  occurred  around  10  m, 
which  is  in  the  vicinity  of  the  Fresnel  radius.  So  the  Fresnel  radius  may  be  considered  as 
a  critical  threshold  for  the  values  of  outer  scale  regarding  the  disturbance  of  a  signal  by 
turbulence.  In  other  words,  since  we  assumed  10  m  as  the  reference  value  for  the 
turbulence  outer  length  scale,  the  correlation  curves  should  begin  to  exhibit  a 
decorrelation  at  this  value.  As  we  increase  this  scale  parameter  to  values  greater  than  the 
Fresnel  radius,  then  the  correlation  curves  may  change  slightly.  But  as  we  decrease  this 
parameter  to  values  less  than  the  Fresnel  radius,  a  dramatic  change  in  coherence  should 
occur,  and  the  correlation  should  increase  considerably. 
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Figures  12  to  14  show  one-dimensional  slices  of  the  turbulent  perturbation  energy 
density  spectrums  in  the  wavenumber  domain  for  different  outer  length  scales.  It  is 
important  to  note  that  around  2  m'\  the  amplitudes  of  the  spectra  are  very  similar. 


Spectral  Density  of  Turbulent  Sound  Speed  Perturbation  (1-D  slice) 


Figure  12.  One-dimensional  energy  density  spectra  for  reference  outer  length  scale  (10 
m)  and  values  below  (2m  and  5m).  Vertical  lines  are  the  Fresnel  radius  values  for  4,  8, 
16,  and  20  kHz,  from  left  to  right,  respectively. 
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Spectral  Density  of  Turbulent  Sound  Speed  Perturbation  (1-D  slice) 


Figure  13.  One-dimensional  energy  density  speetra  for  reference  outer  length  scale  (10 
m)  and  values  above  (25m  and  50m).  Vertical  lines  are  the  Fresnel  radius  values  for  4,  8, 
16,  and  20  kHz,  from  left  to  right,  respectively. 


Spectral  Density  of  Turbulent  Sound  Speed  Perturbation  (1  -D  slice) 


Figure  14.  One-dimensional  energy  density  spectra  for  reference  outer  length  scale  (10 
m)  and  values  above  (25m  and  50m)  on  a  smaller  scale.  Vertical  lines  are  the  Fresnel 
radius  values  for  4,  8,  16,  and  20  kHz,  from  left  to  right,  respectively. 
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Another  important  point  to  make  about  the  influence  of  the  turbulent  outer  length 
scale  is  that  changing  this  value  changes  the  RMS  sound  speed  perturbation  value.  For 
the  reference  turbulent  strength  scaling,  sample  RMS  values  for  the  five  different  outer 
length  cut-off  scales  are  shown  in  Table  4. 


Case 

Ref-2 

Ref-1 

Reference 

Ref+1 

Ref+2 

RMS 

0.217 

0.338 

0.448 

0.627 

0.799 

Table  4.  RMS  Perturbation  of  Sound  Speed  (m/s)  for  Different  Cases  of  Turbulence 
Outer  Length  Scale  Parameter 

3.  Sound  Speed  Profiles 

During  the  ASIAEX,  54  CTD  casts  were  made  by  the  Chinese  research  vessel 
Shiyan-3  from  2  to  6  June,  2001.  The  reference  SSP  is  taken  as  the  average  value  of  these 
measurements.  In  addition  to  the  reference  sound  speed  profile,  three  other  individual 
sound  speed  profiles  were  examined  in  this  analysis.  One  was  chosen  as  being  reasonably 
close  to  the  mean  (reference)  profile,  one  was  chosen  as  being  the  most  deviated  from 
mean  with  a  deeper  thermocline,  and  the  third  was  chosen  as  being  the  most  deviated 
from  mean  with  the  uppermost  thermocline.  These  three  profiles  are  referred  to  as 
Closest  to  Mean  (CTM),  Most  Deviated  Down  (MDD),  and  Most  Deviated  Up  (MDU), 
respectively.  Figure  1 1  shows  how  the  four  sound  speed  profiles  compare. 

Analysis  of  the  influence  of  background  sound  speed  profile  is  more  complex 
than  the  other  environmental  parameters  because  the  variability  introduced  in  the  signal 
coherence  is  not  easily  predictable.  Indeed,  we  initially  assumed  that  the  background 
profile  would  not  significantly  affect  such  short  range  coherence  results.  The  analysis 
presented  in  the  following  chapter,  however,  does  indicate  a  significant  dependence  on 
the  background  profile  and  the  propagation  path  geometry. 

Previous  work  by  Duda  and  Bowlin  (1994)  examined  the  influence  of  sound 
speed  relative  curvature  on  signal  stability.  However,  this  approach  was  examined  in  the 
context  of  long-range,  deep-ocean  environments  and  assumed  to  hold  for  multi-cycle  ray 
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propagation.  Therefore,  it  is  unclear  if  a  similar  measure  of  the  background  sound  speed 
profiles  in  this  environment  will  relate  in  any  way  to  signal  coherence  along  specific  ray 
path  geometries. 


Sound  Speed  Profiles 


Figure  1 1 .  Sound  Speed  Profiles  Used  in  the  Modeling  Analysis 
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IV.  POST  PROCESSING  AND  ANALYSIS 


As  explained  in  Chapter  II,  the  vertieal  spatial  eorrelation  is  used  to  examine  the 
signal  degradation.  If  signal  degradation  is  not  present,  then  the  eorrelation  is  going  to  be 
very  elose  to  unity.  A  perfeet  eorrelation  value  of  1  is  the  expeeted  result  of  an 
autoeorrelation  (the  eorrelation  of  the  signal  with  itself).  All  other  eombinations  for 
eorrelation  are  normalized  aeeording  to  the  autoeorrelation  result,  and  thus  they  are  less 
than  or  equal  to  1 . 

The  eorrelations  are  ealeulated  as  a  funetion  of  the  “depth 
differenee/wavelength.”  The  eenter  frequeney  of  the  signal  and  Cg  =  1500m/ s  are  used 
for  the  wavelength  estimation. 

For  eaeh  parameter,  the  referenee  ease  is  shown  first  followed  by  the  various 
results  for  other  values  of  that  parameter.  Therefore,  some  figures  are  ineluded  several 
times  to  allow  for  eonvenient  eomparison.  “Referenee  Case”  means  that  all  parameters 
are  adjusted  to  their  referenee  values. 

During  the  analysis  of  the  turbulenee  outer  length  seale,  sound  speed  profde,  and 
souree  depth  parameters,  the  referenee  turbulenee  strength  didn’t  ereate  enough 
variability  for  4  kHz  and  8  kHz.  Therefore,  during  the  analysis  of  these  parameters,  the 
referenee  seale  for  these  frequeneies  is  ehanged.  Speeifieally,  the  Ref+2  turbulenee 
strength  parameter  is  used  for  the  4  kHz  analysis,  and  the  Ref+1  turbulenee  strength 
parameter  is  used  for  8  kHz  analysis.  All  graphs  exhibit  four  different  eorrelation  plots 
for  the  different  reeeiver  sub-arrays,  with  the  uppermost  plot  eorresponding  to  the 
shallowest  sub-array  and  the  lowermost  eorresponding  to  the  deepest. 

A,  TURBULENCE  STRENGTH 

The  turbulenee  strength  is  a  parameter  that  direetly  affeets  the  scale  of  turbulence 
in  the  environment.  As  expected,  the  influence  of  this  parameter  is  found  to  depend  upon 
frequency  of  the  propagation.  Examining  Figures  15  through  18,  it  can  be  easily  seen  that 
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as  the  frequency  is  increased,  the  same  turbulence  strength  produces  significantly 
more  variability  in  the  propagation.  This  effect  is  most  noticeable  when  the  frequency 
is  doubled  (from  4  kHz  to  8  kHz,  then  8  kHz  to  16  kHz),  but  is  not  as  evident  above  16 
kHz  in  this  analysis.  Figures  17  and  18  show  that  20  kHz  propagation  has  almost  the 
same  vertical  coherence  as  16  kHz. 

Another  important  result  to  note  is  that  the  coherence  for  the  deepest  sub-array  is 
a  lot  higher  than  the  other  sub-arrays  for  almost  all  cases  of  varying  frequency  or  strength 
parameter.  This  may  be  related  to  the  background  sound  speed  profile  with  which  the 
acoustic  path  interacts.  This  will  be  examined  further  in  a  subsequent  section  using  a 
simple  ray  tracing  model. 

Also  note  that  these  figures  indicate  little  effect  of  the  turbulence  at  4  kHz  when 
the  reference  strength  parameter  is  employed.  As  previously  noted,  all  subsequent 
analysis  at  4  kHz  will  be  performed  using  the  highest  (Ref+2)  turbulence  strength 
parameter.  Similarly  at  8  kHz,  all  subsequent  analysis  will  be  performed  using  the 
second  highest  (Ref+1)  turbulence  strength  parameter. 
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(b)  (c) 


(d)  (e) 

Figure  15.  Effect  of  Turbulence  Strength  Parameter  at  4  kHz. 

(a)  Turbulence  Strength  =  Reference,  (b)  Turbulence  Strength  =  Ref-2,  (c)  Turbulence 
Strength  =  Ref-1,  (d)  Turbulence  Strength  =  Ref+1,  (e)  Turbulence  Strength  =  Ref+2 
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(b) 


(c) 


(d)  (e) 

Figure  16.  Effect  of  Turbulence  Strength  Parameter  at  8  kHz. 

(a)  Turbulence  Strength  =  Reference,  (b)  Turbulence  Strength  =  Ref-2,  (c)  Turbulence 
Strength  =  Ref-1,  (d)  Turbulence  Strength  =  Ref+1,  (e)  Turbulence  Strength  =  Ref+2 
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(b) 


(c) 


(d)  (e) 

Figure  17.  Effect  of  Turbulence  Strength  Parameter  at  16  kHz. 

(a)  Turbulence  Strength  =  Reference,  (b)  Turbulence  Strength  =  Ref-2,  (c)  Turbulence 
Strength  =  Ref-1,  (d)  Turbulence  Strength  =  Ref+1,  (e)  Turbulence  Strength  =  Ref+2 
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(a) 


(b) 


(c) 


(d)  (e) 

Figure  18.  Effect  of  Turbulence  Strength  Parameter  at  20  kHz. 

(a)  Turbulence  Strength  =  Reference,  (b)  Turbulence  Strength  =  Ref-2,  (c)  Turbulence 
Strength  =  Ref-1,  (d)  Turbulence  Strength  =  Ref+1,  (e)  Turbulence  Strength  =  Ref+2 
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B,  TURBULENCE  OUTER  LENGTH  SCALE 

The  analysis  of  the  turbulence  outer  length  scale  confirmed  that  acoustic 
propagation  through  a  turbulent  environment  is  most  sensitive  to  structures  on  the  scale 
of  the  Fresnel  radius.  If  the  range  is  constant,  the  Fresnel  radius  is  a  function  of  frequency 
only.  In  this  model,  the  range  is  kept  constant  and  the  frequency  range  of  4-20  kHz  is 
used.  As  a  result  the  Fresnel  radius  values  vary  from  13.7  m  down  to  6.1  m.  Rather  than 
changing  the  cutoff  scale  for  each  frequency,  kt  =  0.628  1/m  (outer  length  scale  of  /l(=  10 
m)  is  used  as  a  reference  value,  consistent  with  previous  work  (Duda  et  ah,  1998).  Note 
that  this  falls  within  the  range  of  Fresnel  scales  defined.  Thus,  the  influence  of  the  outer 
length  scale  is  considered  here. 

Figures  19  to  22  show  the  effect  of  changes  in  cutoff  scale  at  different 
frequencies.  The  vertical  correlation  is  found  to  exhibit  the  greatest  degradation  (most 
significant  scattering)  for  the  largest  outer  length  scale  value  of  50  m.  This  result  is 
verified  for  all  frequencies. 

A  turbulence  outer  length  scale  of  5  m  barely  produced  any  visible  signal 
degradation,  and  a  cutoff  scale  of  2  m  produced  no  observable  signal  degradation  for  all 
frequencies.  For  these  cutoff  scale  values,  the  turbulence  spectrum  does  not  generate 
enough  perturbation  at  the  Fresnel  scale  for  the  frequencies  considered.  However,  as  the 
cutoff  scale  is  increased  to  10  m,  the  part  of  the  spectrum  generating  perturbations  at  the 
Fresnel  scale  starts  to  increase,  and  the  signal  correlation  tends  to  decrease.  This  result  is 
then  consistent  with  the  expectation  of  the  dominant  effect  of  scattering  at  the  Fresnel 
scale. 

For  4  kHz,  the  Fresnel  scale  is  13.7  m.  When  the  reference  cutoff  scale  of  10  m  is 
applied,  the  perturbations  are  below  the  Fresnel  scale.  Therefore,  a  cutoff  scale  value  of 
25  m  has  a  remarkable  effect  for  this  frequency.  The  same  transition  may  also  be 
observed  for  8  kHz,  since  the  Fresnel  scale  value  for  this  frequency  is  9.7  m,  which  is 
slightly  below  the  reference  cutoff  scale.  The  Fresnel  scale  values  for  16  and  20  kHz  are 
6.8  and  6.1  m,  respectively.  When  the  reference  cutoff  scale  is  applied,  the  perturbations 
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are  already  beyond  the  eritieal  threshold.  Therefore,  inereasing  the  eutoff  seale  to  25  m 
for  these  frequeneies  does  not  eause  a  remarkable  degradation  in  signal  eoherenee. 

Comparison  of  the  analysis  of  8  kHz  versus  16  and  20  kHz  verifies  that  the 
dissipation  of  this  effeet  does  not  oeeur  right  beyond  the  threshold,  but  several  meters 
beyond  the  threshold.  When  a  10  m  eut-off  seale  was  applied  for  8  kHz,  the  perturbations 
are  only  0.3  m  beyond  the  threshold.  However,  these  same  perturbations  are  3.2  and  3.9 
m  beyond  the  threshold  for  16  and  20  kHz,  respeetively.  The  impaet  of  this  ehange  of 
seale  on  the  level  of  perturbation  at  the  Fresnel  seale  is  readily  observed  in  the  speetral 
plots,  displayed  previously  in  Figures  13  and  14.  Though  an  exaet  value  was  not 
ealeulated,  we  eoneluded  that  if  the  perturbations  are  several  meters  greater  than  the 
Fresnel  seale,  further  inerease  of  the  eut-off  seale  doesn’t  signifieantly  affeet  the  signal 
eoherenee.  Figures  19  to  22  elearly  display  the  verifieation  of  these  results. 

Although  larger  eutoff  seales  eould  have  been  used  for  further  analysis,  the 
referenee  value  of  10  m  was  ehosen  for  the  remainder  of  this  thesis.  The  deeision  was 
made  so  that  the  eutoff  seale  was  eonsistent  with  previous  researeh,  and  beeause  this 
seale  did  produee  notieeable  signal  degradation.  Therefore,  the  infiuenee  of  other 
parameters  eould  be  observed  at  this  eutoff  seale. 
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(a) 


(b)  (c) 


(d)  (e) 


Figure  19.  Effect  of  Turbulence  Outer  Length  Scale  at  4  kHz. 

(a)  Turbulence  Scale  =  10  m  (Reference  Case),  (b)  Turbulence  Scale  =  2  m,  (c) 
Turbulence  Scale  =  5  m,  (d)  Turbulence  Scale  =  25  m,  (e)  Turbulence  Scale  =  50  m 
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(a) 


(b) 


(c) 


(d)  (e) 

Figure  20.  Effect  of  Turbulence  Outer  Length  Scale  at  8  kHz. 

(a)  Turbulence  Scale  =  10  m  (Reference  Case),  (b)  Turbulence  Scale  =  2  m,  (c) 
Turbulence  Scale  =  5  m,  (d)  Turbulence  Scale  =  25  m,  (e)  Turbulence  Scale  =  50  m 
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(b) 


(c) 


(d)  (e) 

Figure  21.  Effect  of  Turbulence  Outer  Length  Scale  at  16  kHz. 

(a)  Turbulence  Scale  =  10  m  (Reference  Case),  (b)  Turbulence  Scale  =  2  m,  (c) 
Turbulence  Scale  =  5  m,  (d)  Turbulence  Scale  =  25  m,  (e)  Turbulence  Scale  =  50  m 
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Figure  22.  Effect  of  Turbulence  Outer  Length  Scale  at  20  kHz. 

(a)  Turbulence  Scale  =  10  m  (Reference  Case),  (b)  Turbulence  Scale  =  2  m,  (c) 
Turbulence  Scale  =  5  m,  (d)  Turbulence  Scale  =  25  m,  (e)  Turbulence  Scale  =  50  m 
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c. 


CURVATURE  ANALYSIS 


Previous  work  by  Duda  and  Bowlin  (1994)  examined  the  stability  of  long-range 
aeoustie  ray  propagation  through  internal  wave  fluetuations  in  deep  water.  They  found 
that  the  stability  of  the  rays  propagating  through  various  depths  of  the  waveguide  were 
related  to  mean  values  of  log|f/| ,  where 


U  =  c 


(93) 


and  c  =  c(z)  is  the  depth-dependent  sound  speed  profde.  More  reeent  work  by  Brown,  et 
al.  (2004,  in  press)  also  showed  theoretieally  how  long-range  ray  stability  may  be 
dominated  by  features  of  the  baekground  sound  speed  profde. 

While  both  of  the  aforementioned  analyses  were  eoneemed  with  long-range,  deep 
water  propagation,  the  analysis  presented  here  suggests  a  potentially  similar  influenee  of 
the  baekground  sound  speed  profile.  Therefore,  we  shall  now  examine  a  similar  measure 
of  the  direet  path  propagation  in  the  shallow  water  environment  presented  here. 
Speeifieally,  the  depth-dependent  parameter  U ,  defined  above,  will  be  eomputed  for  the 
four  profiles  examined.  A  simple  ray-traee  model  will  then  be  employed  to  eompute  ray 
trajeetories  through  the  unperturbed  (no  turbulenee)  baekground  profiles  and  the 
subsequent  average  of  /og|f/|  along  the  ray  path  will  be  determined.  While  Duda  and 
Bowlin  (1994)  ineluded  their  internal  wave  perturbations  in  the  ealeulations  of  U,  the 
turbulenee  fluetuations  are  not  expeeted  to  signifieantly  affeet  the  average  of  /og|f/| 
along  the  direet-path  ray  trajeetories. 

In  order  to  eompute  U,  the  diseretely  sampled  sound  speed  profiles  provided  by 
the  measured  data  sets  must  be  fitted  by  some  twiee  differentiable  funetion.  Cubie 
splines  are  employed  here.  Although  Duda  and  Bowlin  (1994)  refer  to  possible  spurious 
values  from  sueh  eubie  polynomials,  this  tends  to  be  a  more  signifieant  problem  when  the 
spaeing  between  measured  sound  speed  values  is  not  uniform  in  depth.  In  the  ease  of  the 
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ASIAEX  data  sets  used  here,  the  depth  sampling  is  uniform,  so  the  use  of  cubic 
polynomials  is  expected  to  be  more  stable. 

As  an  example  of  the  data  fit  achieved  by  the  cubic  spline  technique.  Figure  23 
displays  the  CTM  profile  data  (as  discrete  points)  compared  to  data  interpolated  from  the 
cubic  polynomial.  This  comparison  shows  the  density  of  measured  data  points  relative  to 
the  observed  curvature  in  the  profile,  and  justifies  the  use  of  the  cubic  splines  to  fit  the 
data.  A  direct  calculation  of  log|l7|  is  then  possible.  Figure  24  provides  the  results  of 
these  calculations  for  all  four  sound  speed  profiles  considered. 


Figure  23.  Measured  Sound  speed  Data  Compared  to  Cubic  Spline  Fit  for  CTM  Profile. 
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Log  |U|  (Reference) 


f  10  If 

Log  |U|  (CTM) 


Log  |U|  (MOD) 


Log  |U|  (MDU) 


Figure  24.  Profiles  of  log|f/|  for  the  four  environmental  profiles  examined;  reference 
(upper  left);  CTM  (upper  right);  MDD  (lower  left);  MDU  (lower  right). 


To  compute  ray  trajectories,  a  simple  4*-order  Runge-Kutta  algorithm  was 
employed.  Such  an  algorithm  requires  a  continuous  sampling  of  both  the  sound  speed 
profile  and  its  first  derivative,  both  of  which  are  provided  by  the  cubic  spline 
interpolation.  Although  the  code  can  accurately  handle  smooth  boundary  reflections,  such 
ray  paths  were  not  of  interest  here,  so  only  those  rays  without  boundary  reflections  were 
saved  for  analysis. 

At  each  range  step  along  each  ray  path,  the  value  of  log|f/|  was  then  evaluated. 
By  simply  summing  all  such  values  and  dividing  by  the  number  of  range  steps  in  the 
calculation  out  to  500m,  an  along-ray  average  of  loglf/l  was  computed.  In  Figure  (4), 
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sample  plots  of  ray  paths  in  the  reference  environment  for  three  different  source  depths 
(30m,  50m,  and  70m)  are  displayed.  It  is  important  to  note  that,  in  order  to  compute  a 
well-sampled  profile  of  the  along-ray  mean  of  log|f/|,  a  dense  sampling  of  ray  arrival 

depths  at  the  final  range  (500m)  is  needed.  Because  the  structure  of  the  background 
sound  speed  profiles  tended  to  cause  a  divergence  of  ray  paths  near  the  center  of  the 
water  column  (refer  to  50m  source  depth  ray  paths  in  Figure  25),  various  techniques  were 
used  to  increase  the  density  of  ray  launch  angles  near  0°,  and  the  minimum  and  maximum 
values  of  launch  angles  were  varied  to  insure  arrival  depths  between  the  depths  of  interest 
(about  20m-80m).  Results  of  this  analysis  are  also  displayed  in  Figure  25  for  the  same 
environment  and  source  depths,  but  using  an  altered  code  that  increases  the  density  of  ray 
paths  for  adequate  sampling  of  the  final  profile  of  the  along-ray  mean  of  log|l7| . 
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Figure  25.  Sample  ray  paths  for  the  reference  environment  and  the  corresponding  profiles 
of  the  along-ray  mean  of  log|f/|  for  three  source  depths:  30m  (upper  plots);  50m  (middle 
plots);  70m  (lower  plots). 

The  profiles  of  the  along-ray  mean  of  log|l7|  were  computed  for  all  background 

sound  speed  profde  environments  and  source  depths  considered.  Once  these  profdes 
were  determined,  a  final  average  value  over  the  depth  of  the  sub-arrays  was  then 
computed.  These  unique  values  are  tabulated  and  presented  in  the  following  sections  to 
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compare  with  the  observed  dependence  of  signal  correlation  on  propagation  paths 
through  the  various  environments. 

D,  SOURCE  DEPTH 

The  depth  of  the  source  will  have  a  direet  influence  on  the  propagation  paths 
between  the  source  and  receiver  elements  as  dietated  by  the  sound  speed  profde.  The 
reference  SSP  is  used  for  the  analysis  of  source  depth  influence.  The  reference  SSP  is 
approximately  isospeed  over  the  depth  intervals  0-20  m  and  80-100  m,  and  as  a 
eonsequence  doesn’t  exhibit  a  remarkable  variability  at  these  depths.  There  is  a 
thermocline  between  the  depths  20-80  m,  and  so  propagation  through  this  region 
undergoes  significant  refraction. 

As  the  source  depth  was  increased  from  30  to  40  and  then  50  m,  an  observable 
pattern  existed  in  the  change  of  signal  eoherence  for  all  frequencies.  When  the  source 
was  placed  at  30  m,  the  signal  mostly  propagated  in  the  upper,  isospeed  part  of  the 
waveguide.  These  paths  appear  coincident  with  little  signal  degradation.  As  the  depth  of 
the  source  is  lowered  to  40  m,  the  signal  partly  propagates  through  the  thermocline.  This 
caused  a  slight  decrease  in  the  correlation  for  all  depths  and  frequencies.  When  the 
source  depth  was  changed  to  50  m,  the  sub-arrays  at  45  and  55  m  received  signals  whieh 
traversed  along  the  thermocline.  The  signals  recorded  at  these  sub-arrays  exhibit  a 
remarkable  drop  in  the  signal  coherence.  Verification  of  these  results  is  shown  in  Figures 
26  to  29  (parts  a,  b  and  c  of  figures). 

When  the  source  was  lowered  to  60  m,  signals  at  the  upper  sub-arrays,  35  and  45 
m,  observed  significantly  more  variability  compared  to  the  lower  sub-arrays,  consistent 
with  higher  variability  along  those  paths  propagating  through  regions  of  higher  sound 
speed  gradient.  Placing  the  source  at  a  depth  of  70  m  didn’t  eause  a  noteworthy 
degradation  of  correlation  at  any  receiver  depths  for  all  frequencies  tested  in  this 
research. 

Due  to  this  apparent  correlation  between  background  sound  speed  gradient, 
propagation  paths,  and  signal  coherence,  the  analysis  of  the  variable  source  depth  effect 
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includes  the  effeet  of  SSP  eurvature,  as  described  in  the  previous  seetion.  The  effect  of 
SSP  curvature  was  not  consistent  for  all  results  of  souree  depths.  However,  eonsidering 
that  this  effeet  was  originally  derived  for  deep  and  long-range  propagation,  the 
exeeptions  may  be  due  to  the  shallow  water,  short-range  nature  of  the  problem. 

The  sub-array  averages  of  eorrelation  values  are  calculated  to  provide  a 
eomparison  between  the  signal  coherenee  and  SSP  curvature.  All  of  the  eorrelation 
eurves  generally  exhibited  a  deerease  in  eoherenee  over  the  first  few  wavelengths,  and  a 
slightly  inereasing  coherence  pattern  near  the  last  wavelength.  The  stabilized  region  of 
the  correlation  curve  is  extracted  after  removing  the  correlation  values  corresponding  to 
the  first  five  and  the  last  wavelengths.  These  stabilized  correlation  eurve  values  were  then 
used  to  generate  an  average  eoherenee  value  for  each  sub-array. 

When  there  was  a  noticeable  relative  differenee  between  the  eorrelation  values  of 
receivers,  the  relative  ehange  of  the  SSP  eurvature  parameter  was  eonsistent.  For  the 
souree  depths  of  30  and  40  m,  a  significant  decorrelation  was  not  noted  at  any  depths  and 
frequeneies.  Therefore,  the  infiuenee  of  SSP  curvature  isn’t  analyzed  for  these  souree 
depths.  The  sub-array  averages  of  along-ray  means  of  log|f/|  for  the  reference  sound 
speed  profile  and  depths  of  50,  60,  and  70  m  are  exhibited  in  Table  5.  The  sub-array 
averages  of  eorrelation  values  for  souree  depths  of  50,  60,  70  m  are  exhibited  in  Tables  6 
to  8. 


Receiver 

50  m 

60  m 

70  m 

35-40  m 

3.54 

3.33 

3.26 

45-50  m 

3.43 

3.13 

2.97 

55-60  m 

3.06 

2.82 

3.00 

65-70  m 

2.76 

3.00 

3.03 

Table  5.  Sub-Array  Averages  of  Along-Ray  Means  of  log|f/|. 
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Receiver''''*^^ 

4  kHz 

8  kHz 

16  kHz 

20  kHz 

35-40  m 

0.91 

0.88 

0.94 

0.92 

45-50  m 

0.75 

0.82 

0.85 

0.89 

55-60  m 

0.85 

0.88 

0.89 

0.88 

65-70  m 

0.95 

0.97 

0.98 

0.98 

Table  6.  Sub-Array  Averages  of  Correlation  Values  for  Souree  Depth  =  50  meters. 


Receiver''''*''''''^^ 

4  kHz 

8  kHz 

16  kHz 

20  kHz 

35-40  m 

0.86 

0.86 

0.90 

0.85 

45-50  m 

0.78 

0.76 

0.84 

0.89 

55-60  m 

0.94 

0.94 

0.97 

0.93 

65-70  m 

0.94 

0.95 

0.97 

0.97 

Table  7.  Sub-Array  Averages  of  Correlation  Values  for  Souree  Depth  =  60  meters. 


Receiver''''**'''*^^ 

4  kHz 

8  kHz 

16  kHz 

20  kHz 

35-40  m 

0.84 

0.89 

0.92 

0.87 

45-50  m 

0.87 

0.89 

0.89 

0.93 

55-60  m 

0.92 

0.86 

0.93 

0.89 

65-70  m 

0.95 

0.93 

0.93 

0.93 

Table  8.  Sub-Array  Averages  of  Correlation  Values  for  Source  Depth  =  70  meters. 


As  the  average  along-ray  mean  values  of  log|f/|  corresponding  to  each  sub-array 
increase,  the  signal  coherence  is  expected  to  decrease.  This  negative  effect  of  SSP 
curvature  may  be  observed  by  cross-examining  Tables  5  to  8.  There  are  exceptions  of  this 
effect,  which  again  may  be  due  in  part  to  the  shallow  water,  short-range  nature  of  the 
problem. 
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(a) 


(b)  (c) 


(d)  (e) 


Figure  26.  Effect  of  Source  Depth  Variation  at  4  kHz. 

(a)  Source  Depth  =  50  m  (Reference  Case),  (b)  Source  Depth  =  30  m,  (c)  Source  Depth  = 
40  m,  (d)  Source  Depth  =  60  m,  (e)  Source  Depth  =  70  m 
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(b)  (c) 


(d) 


(e) 


Figure  27.  Effect  of  Source  Depth  Variation  at  8  kHz. 

(a)  Source  Depth  =  50  m  (Reference  Case),  (b)  Source  Depth  =  30  m,  (c)  Source  Depth  = 
40  m,  (d)  Source  Depth  =  60  m,  (e)  Source  Depth  =  70  m 


60 


(a) 


(b)  (c) 


(d)  (e) 

Figure  28.  Effect  of  Source  Depth  Variation  at  16  kHz. 

(a)  Source  Depth  =  50  m  (Reference  Case),  (b)  Source  Depth  =  30  m,  (c)  Source  Depth  = 
40  m,  (d)  Source  Depth  =  60  m,  (e)  Source  Depth  =  70  m 
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Figure  29.  Effect  of  Source  Depth  Variation  at  20  kHz. 

(a)  Source  Depth  =  50  m  (Reference  Case),  (b)  Source  Depth  =  30  m,  (c)  Source  Depth  = 
40  m,  (d)  Source  Depth  =  60  m,  (e)  Source  Depth  =  70  m 
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E.  SOUND-SPEED  PROFILE 

Four  different  baekground  sound  speed  profiles  were  examined  in  this  researeh. 
Details  of  these  profiles  were  deseribed  in  Chapter  III.  During  the  following  analysis  of 
the  effeet  of  baekground  SSP,  the  souree  depth  was  fixed  at  50  m. 

As  noted  previously,  the  referenee  SSP  is  nearly  isospeed  in  the  intervals  0-20  m 
and  80-100  m,  and  eontains  a  thermoeline  between  the  depths  20-80  m.  The  signal 
eoherenee  was  already  observed  to  suffer  the  most  degradation  when  the  signal 
propagates  through  the  thermoeline.  Sub-arrays  loeated  at  45-50  and  55-60  m  were  found 
to  reeeive  the  most  degraded  signals  when  the  referenee  SSP  is  used.  This  observation  is 
verified  for  all  frequeneies  and  the  results  are  exhibited  in  Figures  30  to  33  (upper  left 
parts  of  all  figures). 

The  CTM  sound  speed  profile  has  depth  properties  similar  to  the  referenee  SSP. 
However,  the  CTM  profile  has  more  variability  in  the  thermoeline  region.  This  inerease 
of  variability  in  sound  speed  is  found  to  produee  even  more  degradation  in  signal 
eoherenee.  A  deerease  in  eorrelation  at  all  depths  is  observed.  It  must  be  noted  that  the 
signals  on  the  sub-arrays  loeated  at  45-50  and  55-60  m  still  exhibit  the  greatest 
degradation  eompared  to  other  sub-arrays,  presumably  sinee  they  are  loeated  in  the 
middle  of  thermoeline  region.  These  results  may  be  eonfirmed  by  looking  at  Figures  30 
to  33  (upper  right  parts  of  all  figures). 

The  MDD  sound  speed  profile  is  mostly  isospeed  over  the  depths  0-65  m  (with 
the  exeeption  of  a  single  “jump”  around  30  m).  Therefore,  based  on  the  previous 
findings,  we  may  expeet  that  signals  from  a  50  m  souree  will  not  be  as  severely  degraded 
as  they  are  with  the  referenee  profile,  sinee  most  of  the  propagation  oeeurs  in  this  upper, 
isospeed  portion  of  the  profile.  The  analysis  of  the  results  using  the  MDD  sound  speed 
profile  eonsistently  verified  this  assumption  and  the  results  are  displayed  in  Figures  30  to 
33  (lower  left  parts  of  all  figures). 

The  MDU  sound  speed  profile  has  effeetively  no  upper  isospeed  region.  The 
thermoeline  exists  from  around  5  m  down  to  around  65  m,  below  whieh  the  profile  is 
nearly  isospeed.  The  overall  variability  of  the  MDU  sound  speed  profile  appears  less  than 
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the  CTM  sound  speed  profile.  Therefore,  it  may  be  antieipated  that  the  results  of  these 
two  sound  speed  profdes  from  a  50  m  souree  would  be  similar,  with  the  MDU  profde 
having  a  weaker  effeet  on  signal  eoherenee.  The  model  verified  these  assumptions  and 
the  results  are  shown  in  Figures  30  to  33  (upper  right  parts  of  all  figures). 

Similar  to  analysis  of  souree  depth,  effeet  of  SSP  eurvature  is  also  analyzed  for 
ehanges  in  the  baekground  sound  speed  profile.  The  most  eonsistent  results  were  noted 
for  the  referenee  SSP  and  the  CTM  sound  speed  profile.  The  sub-array  averages  of  along- 
ray  means  of  log|f/|  for  the  50  m  souree  depth  in  the  various  sound  speed  profiles  are 
tabulated  in  Table  9.  The  sub-array  averages  of  eorrelation  values  for  the  various  sound 
speed  profiles  are  displayed  in  Tables  10-11. 


Receiver'''''*''''^^ 

Reference 

CTM 

MDD 

MDU 

35-40  m 

3.54 

4.60 

4.67 

4.62 

45-50  m 

3.43 

4.86 

5.30 

4.27 

55-60  m 

3.06 

4.05 

4.38 

4.26 

65-70  m 

2.76 

3.69 

3.87 

4.26 

Table  9.  Sub-array  averages  of  along-ray  means  of  log|f/|  for  the  source  at  50  m  depth. 


Receiver'''''''''^^ 

4  kHz 

8  kHz 

16  kHz 

20  kHz 

35-40  m 

0.91 

0.88 

0.94 

0.91 

45-50  m 

0.75 

0.82 

0.85 

0.75 

55-60  m 

0.85 

0.88 

0.89 

0.85 

65-70  m 

0.95 

0.97 

0.98 

0.95 

Table  10.  Sub-Array  Averages  of  Correlation  Values  for  the  reference  SSP. 
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Reeeiver'''''''''*''^^ 

4  kHz 

8  kHz 

16  kHz 

20  kHz 

35-40  m 

0.94 

0.94 

0.94 

0.88 

45-50  m 

0.72 

0.71 

0.73 

0.78 

55-60  m 

0.67 

0.65 

0.74 

0.74 

65-70  m 

0.79 

0.85 

0.88 

0.90 

Table  11.  Sub-Array  Averages  of  Correlation  Values  for  the  CTM  sound  speed  profile. 


As  the  average  along-ray  mean  values  of  log|[/|  eorresponding  to  eaeh  sub-array 
inerease,  the  signal  eoherenee  may  be  expeeted  to  deerease.  The  results  of  the  analysis 
presented  here  verify  this  expeetation  with  some  exeeptions.  The  results  may  be  best 
understood  by  eross-examining  Tables  9  to  11.  The  exeeptions  are  possibly  due  to  the 
short-range  nature  of  the  propagation  examined  in  this  work. 
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(a) 


(b) 


(c)  (d) 

Figure  30.  Effect  of  Background  Sound-Speed  Profde  at  4  kHz. 

(a)  SSP  =  Reference  Case,  (b)  SSP  =  CTM,  (c)  SSP  =  MDD,  (d)  SSP  =  MDU 
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(c) 


(d) 


Figure  31.  Effect  of  Background  Sound-Speed  Profde  at  8  kHz. 

(a)  SSP  =  Reference  Case,  (b)  SSP  =  CTM,  (c)  SSP  =  MDD,  (d)  SSP  =  MDU 
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(a) 


(b) 


Figure  32.  Effect  of  Background  Sound-Speed  Profde  at  16  kHz. 

(a)  SSP  =  Reference  Case,  (b)  SSP  =  CTM,  (c)  SSP  =  MDD,  (d)  SSP  =  MDU 
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(C)  (d) 


Figure  33.  Effect  of  Background  Sound-Speed  Profde  at  20  kHz. 

(a)  SSP  =  Reference  Case,  (b)  SSP  =  CTM,  (c)  SSP  =  MDD,  (d)  SSP  =  MDU 
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V.  SUMMARY 


A,  CONCLUSIONS 

The  motivation  of  this  research  was  to  develop  a  better  understanding  of  the 
influence  of  turbulence  and  environmental  variability  on  broadband  acoustic  coherence. 
The  MMPE  model  was  used  to  simulate  the  effects  of  a  turbulent  environment  on  such 
broadband  acoustic  propagation.  Realistic  turbulence  perturbations  based  on  well 
established  statistical  models  were  incorporated  into  the  MMPE  model,  making 
investigations  of  such  effects  on  signal  coherence  easily  accessible  for  a  variety  of 
environmental  factors  and  source/receiver  geometries.  The  background  sound  speed 
profiles  and  the  model  geometry  were  similar  to  a  component  of  one  of  the  East  China 
Sea  experiments  conducted  during  ASIAEX.  Effects  of  turbulence  strength,  turbulence 
cut-off  scale,  background  sound  speed  profde,  frequency  and  source/receiver  geometries 
on  short  range,  direct-path  propagation  were  examined. 

Vertical  correlation  was  used  to  quantify  the  coherence  of  the  received  signal  at 
four  vertical  sub-arrays  located  at  a  fixed  range  from  the  source.  Each  case  was  simulated 
10  times  to  simple  statistics  (mean  and  standard  deviation)  of  the  signal  decorrelation  due 
to  turbulence. 

It  was  observed  that  the  turbulence  strength,  a  parameter  that  directly  affects 
the  RMS  magnitude  of  the  turbulence  perturbation,  exhibited  a  positive  correlation  with 
signal  degradation  for  all  frequencies  up  to  20  kHz,  as  expected.  Essentially,  this  simply 
states  that  the  stronger  the  sound  speed  perturbations  are,  the  stronger  the  acoustic  signal 
degradation  will  be. 

The  analysis  of  the  turbulence  outer  length  scale  also  confirmed  the  expected 
result  that  acoustic  propagation  through  a  turbulent  environment  is  most  sensitive  to 
structures  on  the  scale  of  the  Eresnel  radius.  When  the  perturbations  are  already  greater 
than  the  Eresnel  scale,  further  increase  in  cut-off  scale  (even  with  increases  in  total  RMS 
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perturbation  strength)  didn’t  affect  the  signal  coherence.  An  exact  value  was  not 
calculated  for  the  critical  excess  value  beyond  which  the  coherence  is  not  affected 
remarkably.  However,  it  was  observed  that  if  the  turbulence  cut-off  scale  is  several 
meters  greater  than  the  Fresnel  scale,  then  further  increases  of  the  cut-off  scale  do  not 
produce  a  noteworthy  effect  on  the  signal  coherence. 

Source-receiver  depth  and  background  SSP  were  found,  somewhat  surprisingly, 
to  have  a  significant  influence  on  the  signal  coherence.  This  was  not  expected  because 
there  was  no  direct  relationship  between  the  background  profile  and  the  structure  of  the 
turbulence,  and  because  the  propagation  distances  were  so  short  that  such  influences 
were  not  expected  to  have  enough  time  to  “build  up”  along  the  path.  It  was  confirmed 
that  propagation  through  a  thermocline  has  greater  (negative)  effect  on  signal  coherence 
than  the  isospeed  regions  of  the  SSP.  Thus,  the  signal  stability  appears  to  be  strongly 
dependent  upon  the  propagation  path  through  the  background  profde  (without 
turbulence). 

Due  to  these  observations,  an  attempt  was  made  to  compare  these  results  with 
similar  findings  of  signal  stability  dependence  in  deep  ocean,  long-range  propagation  on 
sound  speed  profde  curvature.  The  analysis  of  the  SSP  curvature  performed  here  was 
partly  consistent  with  the  findings  of  the  signal  coherence.  Most  notably,  when  there 
were  significant  relative  differences  between  the  correlation  values  of  the  sub-arrays,  the 
relative  differences  of  the  mean  curvature  parameter  along  the  ray  paths  were  consistent. 

B,  RECOMMENDATIONS  AND  FUTURE  WORK 

The  findings  of  the  significance  of  the  structure  of  the  background  sound  speed 
profile  were  rather  surprising,  but  of  considerable  significance.  Such  findings  could  result 
in  changes  in  the  operation  of  various  underwater  acoustic  systems.  For  example,  tactical 
systems  (e.g.,  mine-hunting  UUVs)  may  be  able  to  improve  their  detection  capabilities  by 
determining  their  optimal  operating  depths  based  on  in  situ  measurements  of  the 
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background  profile.  Depth  variable  aeoustie  eommunieation  systems  may  also  be  able  to 
optimize  their  performanee  in  a  similar  fashion. 

While  the  analysis  of  SSP  eurvature  performed  in  this  thesis  showed  some 
eonsisteney  with  the  signal  eoherenee  results,  it  is  important  to  reeall  that  this  measure 
was  originally  derived  for  deep  and  long-range  propagation.  It  is  possible,  therefore,  that 
another  measure  of  the  interaetion  of  the  propagation  paths  with  the  baekground  profile 
may  be  better  suited  for  predietion  of  signal  degradation.  Future  work  may  then  examine 
other  such  measures. 

It  should  also  be  noted  that  this  analysis  assumed  a  uniform  turbulenee 
perturbation  field.  The  effeets  of  depth  or  range  variability  of  the  turbulenee  is  unknown, 
as  is  any  obvious  eorrelation  with  depth  or  range  variability  of  the  baekground  profile. 
More  sophistieated  oeeanographie  models  of  turbulenee  would  be  needed  to  generate  the 
neeessary  realizations  of  turbulent  fiuetuations  in  order  to  study  these  effeets. 

Finally,  although  this  modeling  effort  was  based  on  the  experimental  geometry 
employed  during  the  ASIAEX  studies  in  the  East  China  Sea,  the  measured  data  sets  did 
not  inelude  any  direet  measurements  of  turbulenee  struetures.  It  would  be  of  great 
benefit  in  a  future  experiment  to  simultaneously  measure  both  signal  eoherenee  and  water 
eolumn  fine-seale  strueture.  The  validity  of  this  numerieal  approaeh  may  then  be 
eonfirmed,  and  further  numerieal  studies  eould  then  be  performed  with  eonfidenee.  Our 
understanding  of  the  impaet  of  small-seale  oeean  fiuetuations  on  various  aeoustie 
systems’  performanee  eould  then  be  enhaneed,  and  future  system  designs  eould  be 
improved. 
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